!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! calcaulate g, Rg, Cp for any subroutine
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! coded by : B. G.
! created  : 2015-06-17
! revised  :
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! nomenclature:
!   Rsl     gas constant at SL, 287.05......
!   b0~b5   constants
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

subroutine GasCons(T,Cp)
   use VariableDef
   implicit none
   real*8::T,Rsl,Cp,Cv
   real*8::b0,b1,b2,b3,b4,b5
   
   Rsl=287.05
   b0=1.817160e3
   b1=9.890504e2
   b2=-9.595592e-3
   b3=1.041469e-4
   b4=-4.433065e-8
   b5=5.879263e-12
   
   Cp=b1+2*b2*T+3*b3*T**2+4*b4*T**3+5*b5*T**4
   Cv=b1-Rsl+2*b2*T+3*b3*T**2+4*b4*T**3+5*b5*T**4
   g=Cp/Cv
   Rg=Cp/(g/(g-1))

end subroutine GasCons
